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^ ■ Abstract 

CN . A few-body formalism is applied for computation of two different three-charge-particle systems. 



The first system is a collision of a slow antiproton, p, with a positronium atom: Ps= (e + e~) — a 
bound state of an electron and a positron. The second problem is a collision of p with a muonic 
muonium atom, i.e. true muonium — a bound state of two muons one positive and one negative: 
Ps^ = (// + / u~). The total cross section of the following two reactions: p + (e + e~) — > H + e~ and 
p+(^ + ^~) — > H^ + u~ , where H = (pe + ) is antihydrogen and ELj = (p^ + ) is a muonic antihydrogen 
^* atom, i.e. a bound state of p and // + , are computed in the framework of a set of coupled two- 

i ^j , component Faddeev-Hahn-type (FH-type) equations. Unlike the original Faddeev approach the 

£NJ ' FH-type equations are formulated in terms of only two but relevant components: ^i and \l/2, of 

>' 

the system's three-body wave function *&, where ^ = ^i + ^2- In order to solve the FH-type 

m , 

equations \l/i is expanded in terms of the input channel target eigenfunctions, i.e. in this work 

(N ■ 

tj- | in terms of, for example, the (u + u~) atom eigenfunctions. At the same time \&2 is expanded in 

o 

terms of the output channel two-body wave functions, that is in terms of H M atom eigenfunctions. 
Additionally, a convenient total angular momentum projection is performed. This procedure leads 



^ \ to an infinite set of one-dimensional coupled integral-differential equations for unknown expansion 

H ■ 

functions. Since the two-body targets are treated equally and the accurate asymptotes of ^1 and 

^2 are provided, the solution of the FH-type equations avoids the over-completeness problem. 

Results for better known low-energy (jT transfer reactions from one hydrogen isotope to another 

hydrogen isotope in the cycle of muon catalyzed fusion (/iCF) are also computed and presented. 
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I. INTRODUCTION 



To date, non-relativistic quantum-mechanical Coulomb few-body problems have a long 
research history. In fact, the first works dealing with quantum few-charged particle systems 
appeared during the early stages of quantum-mechanics [l| . This is because these types of 
problems pose significant fundamental theoretical and practical importance in nuclear and 
atomic-molecular physics. The few-body Coulomb problem is of considerable importance 
in the cycle of //CF (cold fusion) [2], in cases where a muonic few-body system experi- 
ences a strong interplay between Coulomb and nuclear forces involving heavy nuclei, for 
instance, the (dtyU.) + molecular ion. Further, it would be worth mentioning that there are 
modern antimatter physics problems, that involve few-body systems, such as antihydro- 



gen [3H8| /protonium 9l-|l3J| formation reactions, low energy p+H^" collisions 14(, H+H 2 
quenching 15|, LL6], H+H annihilation reactions 17j], and p+ 4 ^ 3 ^He antiprotonic helium atom 
(atomcule) formation 18|, Il9] just to name a few. It follows from the charge conjugation, 
parity, and time reversal (CPT) symmetry of quantum electrodynamics that a charged par- 
ticle and its antiparticle should have equal/opposite charges, equal masses, lifetimes, and 
gyromagnetic ratios. Moreover, the CPT symmetry predicts that hydrogen and antihydro- 
gen atoms should have identical spectra. New experiments are in progress to test these 
fundamental laws and theories of physics involving antiparticles as well as antimatter in 
general. Therefore, future experimentalists plan to test whether H and H have such prop- 
erties. In such experiments, it would be important to have a certain quantity of H atoms 
at low kinetic energies, ideally at the rest: T ~ K 20t |21|. Using this perspective we de- 



velop a quantum-mechanical approach which would be reliable at low and very low collision 
energies, i.e. when the quantum-mechanical few-body dynamics of three Coulomb particles 
becomes important. The method is formulated for arbitrary masses of the particles, that is 
when the dynamics of lighter and heavier particles are not separated from each other. 



The author of the book 22j pointed out that the muonic antihydrogen atom, H^, could be 
even a better choice to check the CPT law than the usual antihydrogen atom. This is because 
the size of this atom is ~ 207 times smaller than the size of a normal H atom. Therefore, 
as mentioned in 22(: the short range CPT violating interaction with an extremely heavy 
boson can be easily detected within the system. This idea appears be extremely interesting, 
and therefore it would be useful to compute the formation cross sections and rates of H M at 



low energy collisions, for example, from ~1 eV down to ~ 10~ 5 eV. Thus, in this work we 
consider the following three-body reactions of antihydrogen H and muonic antihydrogen H M 
formation: 

p+(e+e-) ls ^H + e-, (1) 

p+(/i + /i")is^H At + /i^. (2) 

At such low energies the quantum-mechanical Coulomb few-body dynamics become im- 
portant, especially in the case of heavy charge transfer, i.e. /i + . Also, it would be quite 
appropriate to mention that exotic atomic and antiatomic systems like a true muonium 
atom, (fi + fi~), or a simple muonic hydrogen atom, H M =(p + /x~), are always of great interest 



in nuclear, atomic and few-body physics [23H25J . For instance, recently the authors of works 

26L l27| | have considered an interesting problem: the production of (yU + /i~). This is the 

smallest pure QED atom with the Bohr radius only ~512 fm. So far ( , ^ + /i~) has never been 



observed. Next, in the recent works 



28 



29| | the proton-radius puzzle 30(] was considered 



from few-body and muonic physics perspectives. Another three-charge-particle reaction of 



H M formation was considered in the works 



22 



25| too: 



p + Mu -»■ H„ + e" 



(3) 



Here, Mu is the muonium atom, i.e. a bound state of a positive muon /i + and an electron: 



Mu=(fi + e-). 



reaction 
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3l| 



his is a very interesting and complex example of a heavy charge transfer 



In nuclear physics, involving applications related to three-body systems the few-body 
Faddeev and Alt-Grassberger-Sandhas (AGS) equations |32|-|35| are frequently employed. 
These equations are equal to the Schrodinger equation, but formulated for the three-body 
wave function components and therefore have the correct physical asymptotes. However, 
in the case of three-charged particle systems the kernels of the original integral Faddeev 
equations in momentum space lose their compactness due to Coulomb long-range interactions 
33| . This limitation has become the most serious obstacle in the practical application of the 
original Faddeev equation to few-body systems with pure Coulomb interactions. Therefore, 
on one hand the Faddeev and AGS equations are the most rigorous attempt to provide 
a basis for detailed few-body numerical computations, but on the other hand they have 
not been used much to date because they have been regarded as too complex to solve when 
used in Coulomb scattering problems. This limitation has led to various alternative methods. 



Among the most popular is a well known method based on the Born-Oppenheimer adiabatic 



model 



and improved adiabatic approximation [37J. The approach has been applied to 
many systems in atomic and //-atomic physics for over many years. Very accurate variational 
calculations have been applied to selected three-body Coulomb systems in muon catalyzed 



fusion cycle [38], and in H formation reactions, see for instance [3[. Here it would also be 
useful to mention important Coulomb few-body calculations based on the adiabatic hyper- 
spherical method p, |39], coordinate-space Faddeev equation approach in three dimensions 



40| and within a hyper-spherical function expansion formalism 41]. There are also newer 



developments in the field we would also like to cite 421444 ] . 

In the current work however we apply a different few-body approach based on a set of cou- 



pled two-component FH-type equation formalism [45144 71] . The next section represents the 
notation pertinent to the three-charged-particle system (123) shown in Fig. 1, the original 
equations, boundary conditions, detailed derivation of the set of coupled one- dimensional 
integral-differential equations suitable for a numerical computation and the numerical com- 
putational approach developed in this work. Sec. Ill includes new results, conclusions, and 
Sec. IV includes Appendix. The atomic units, i.e. e = h = m e = 1, are used in the case of 



the e + transfer reaction, and the muonic units, i.e. e = H 



m. 



1, are used in the case of 



the /i + //i transfer reactions, where m M = 206.769 m e is the mass of the muon. 



II. FEW-BODY TREATMENT 



In the case of three charged particles (123), two positive and one negative, only two 
asymptotic configurations are possible below the breakup threshold. This situation is shown 
in Fig. 2, for instance, for the (p n~ fi + ) system. It suggests to write down a set of two 
coupled equations for Faddeev- type components of the system's wave function 45[. These 
equations are commonly called Faddeev-Hahn-type (FH-type) equations [46j. In this work 
we shall consider a method based on an integral-differential equation approach [47J applied 
to Coulomb three-body systems. To solve these equations, a modified close coupling method 
is applied. This procedure leads to an expansion of the three-body system wave function 
components into eigenf unctions of the subsystem Hamiltonians, providing an infinite set of 
one-dimensional integral-differential equations [47]. Within this formalism the asymptotic 
of the full three-body wave function contains two parts corresponding to two open channels. 



In this work we consider different Coulomb three-body systems with arbitrary masses, i.e. 
the masses of the charged particles are taken as they are. We do not apply any type of 
adiabatic approximations, when the dynamics of heavy and light parts of the system are 
separated. Therefore, this dynamical method works for both a light charge, e + , transfer 
reaction at ultra- low energies and for heavy charge transfer, p + , as well. 

A. FH-type equation approach 

Let us define the system of units to be e = h = m 3 = 1 and denote antiproton p by 1, 
a negative muon p~ by 2, and a positive muon p + by 3. Before the three-body breakup 
threshold two cluster asymptotic configurations are possible in the three-body system, i.e. 
(23) — 1 and (13)— 2 being determined by their own Jacobi coordinates {rj 3 ,pk} as shown in 
Figs. 1, 2: 

r j3 = r 3 - fj, (4) 

Pk = (5 + mjfj)/(l + rrij) - r k , (j ^ k = 1,2). (5) 

Here r*g, m^ are the coordinates and the masses of the particles £ = 1, 2, 3 respectively. This 
suggests a Faddeev formulation which uses only two components. A general procedure to 
derive such formulations is described in work |45(. In this approach the three-body wave 
function is represented as follows: 

|^) = *i(r 23 ,Pi) + *2(ri3,p 2 ), (6) 

where each Faddeev-type component is determined by its own Jacobi coordinates. Moreover, 
^i (^23 , pi) is quadratically integrable over the variable r 23 , and ^2(^13, P2) over the variable 
r*i3. To define \^i), {I = 1,2) a set of two coupled Faddeev-Hahn-type equations can be 
written: 

(e-H - V 23 (f 23 )^ 1 (f 23 ,p 1 ) = (V 23 (f 23 ) + V 12 (f 12 ))* 2 (r 13 ,p 2 ), (7) 

(e-H - V 13 (f 13 )^ 2 (f l3 ,p 2 ) = (V 13 (r 13 ) + K 2 (fi 2 ))*i(f 2 3,pi). (8) 

Here, Hq is the kinetic energy operator of the three-particle system, Vijfrij) are paired in- 
teraction potentials (i ^ j = 1,2,3), E is the total energy. The constructed equations 
satisfy the Schrodinger equation exactly. For the energies below the three-body break-up 



threshold they exhibit the same advantages as the Faddeev equations [33| , because they are 
formulated for the wave function components with the correct physical asymptotes. To solve 
the equations a close-coupling method is applied, which leads to an expansion of the sys- 
tem's wave function components into eigenfunctions of the subsystem (target) Hamiltonians 
providing with a set of one-dimensional integral-differential equations after the partial-wave 
projection. A further advantage of the Faddeev-type method is the fact that the Faddeev- 
components are smoother functions of the coordinates than the total wave function. Also, 
the Faddeev decomposition avoids overcompleteness problems, because two-body subsys- 
tems are treated in an equivalent way, and the correct asymptotes are guaranteed. Next, 



based on Merkuriev 



33 



34J the three-charge-particle scattering wave function and its all 
three components should have the following general asymptotic form: 

00 piknpk 

Mrjs, Pk) ~ e^VxfoOfci + £ Af e %Q Pk )-—<p n (f j3 ) + 

p k ->-+co * — ' Oh 

n 

~ pik'mPi i(^E P +W c (p,E) 

J2AZ(n Pi )^ Vm (f k3 ) + B(n 5 ) m . (9) 

Here, e tklZ ipi(fj3) is the incident wave, y? n (r} 3 ) the n— th bound-state wave function of the 
pair (j3), k n = (E — £ n ) 1//2 , £n is the binding energy of the (j'3), A el l ex (VL Pk ), A tr (VL Pj ) 
and B(£l§) are amplitudes of elastic/inelastic, transfer and breakup channels respectively, 
Pe, = (p, f2 5 ) is the three-body hyperradius and W c (p, E) is the three-body Coulomb phase 



34l | . For lower energy collisions when E < E t h r , where E t hr is the three-body break-up 
threshold, the expression (Q becomes simpler, i.e. without the last term: B(Q 5 ) = 0. 

Therefore, in the current work for low energy collisions each Faddeev type component 
corresponds to only one determined channel. For example, for the elastic and for the charge 
transfer channel we have: 

00 ik p\ 

*i(f 23 ,Pi) ~ e^^j + ^if^)— ^(f 23 ), (10) 

n 
00 ik' rn 

Mn 3 ,P2) ~ E A -(^) ^(ns). (11) 

m 

It is easy to see that the asymptotic behavior of the total wave function ([6]) becomes similar 
to equation (jSJ). 

In addition, we would like to point out, that the few-body FH-type equation approach 
©-© is a quite flexible method. For example, let us briefly consider a muon transfer 
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reaction in the following low energy collision: 

Li 3+ + (p+/Ois^(Li A 2 t + ) ls + P + . (12) 

This reaction has a strong, pure Coulomb interaction in the output channel. This circum- 
stance can be taken into account by adding a distortion potential into the FH equations 
0-dED, i.e. U(p 2 ) = [Z\ — l)Z 2 /p2, where Z x is the charge of Li 3+ and Z 2 {=1) is the charge 
of the hydrogen isotope: 

(E - h 23 (f 23 ) - f l (p l )^ 1 (r 23 ,p 1 ) = (V 23 (f 23 ) + V 12 (f 12 ) - (Zl ~ 1)Z2 )^(f 13 ,p 2 ), (13) 

(E - h 13 (f 13 ) - f 2 (p 2 ) - (Zl - 1)Z2 )^^(r 13 ,p 2 ) = (v 13 (f l3 ) + V 12 (f 12 )^ 1 (f 23 ,p 1 ). (14) 

The two coupled equations satisfy the Schrodinger equation exactly, i.e. when the Eqs. 
(TT3|) - (fl4|) are added to each other the two distortion potential terms vanish. Also, in these 
equations we identify the target hamiltonians h 23 (f 23 ) and hi 3 (fi 3 ), as well as the kinetic 
energy operators, i.e. Ti(pi) and T 2 (p 2 ). One can see, that by converting the differential 
Eq. 014p into an integral equation we can obtain a Coulomb Green function over the Jacobi 
coordinate p 2 in the output channel. The Green function provides the physically correct 
Coulomb asymptotic for the component \^ 2 ). The component |\Pi) carries the asymptotic 
behavior for the elastic and inelastic channels and the component \^ 2 ) carries the Coulomb 
asymptotic behavior in the transfer channel, that is: 

ifc (1) oi 
/-j_\ p t,r "n l J i- 

*i(r 23 ,Pi) ~ e^ VifoO + J2 A n' m ( n pi) <Pn(?23), (15) 

n ' 

__^ p i(k$ p a -irl/2+ri-r)/2kg ) In2fc^'p2) 

* 2 C (ri3,p 2 ) ~ E A ^(^) Vmins). (16) 



P2^r+CO * ' p 2 

ml 



1.W 



Here, e lkl z ^\{r 23 ) is the incident wave, ip n (fj 3 ) the n-th excited bound-state wave function of 



the pair (j3), kn = \/2M i (E — E n 3 ), with Mf l = m^ 1 + (1 + rrij)' 1 , where m ; represents 
the masses of the heavy particles p and Li 3+ . Here En is the binding energy of (j3), 
i ^ j — 1,2, v4 cl / in (f2 Pl ) and A tr (Q P2 ) are the scattering amplitudes in the elastic/inelastic 
and transfer channels. The Coulomb parameters in the second transfer channel are: 77 = 
argr(/ + 1 + ir] /2km) and r\ = 2M 2 \Z\ — l)/k n ■ One can see, that this approach simplifies 
the solution procedure and provides the correct asymptotic behavior for the solution below 
the three-body breakup threshold. Further, the few-body method has been successfully 



applied to different three-body muon transfer reactions 46| 



B. Obtaining an infinite set of coupled integral-differential FH-type equations 

Now, let us present the equations flZ])-© in terms of the adopted notation 

(E + ^f k A P, + o^ A ^3 - Vf^Vtf&Pk) = {V j3 + V jk )^ v {r k3 ,p^ (17) 

here i ^ %' — 1, 2, M^ 1 = m^ 1 + (1 + m J )" 1 and /17 1 = 1 + m 7 . In order to separate angular 
variables, the wave function components \l/j are expanded over bipolar harmonics: 

{Y x (p) ® >^) W = ^CSL^VOSJlWCr). (18) 



fj.ru 



where p and r are angular coordinates of vectors p and r; C^( m are Clebsh- Gordon coeffi- 



cients; Y/ m are spherical functions |48j. The configuration triangle of the particles (123) is 
presented on the Fig. 1 together with the Jacobi coordinates {r 2 3,pi} and {ri3,p 2 } and an- 
gles between them. The centre-off-mass of the whole three-body system is designated as O. 
The centre-off-masses of the two-body subsystems (23) and (13) are 0\ and O2 respectively. 
Substituting the following expansion: 

*i(*V3, Pk) = J2 *LMXi(Pk, r jS ) {Yx(p k ) <8> Yf(f j3 )} LM (19) 

LMXl 

into ( 1T7J) . multiplying this by the appropriate biharmonic functions and integrating over the 
corresponding angular coordinates of the vectors ?j 3 and p k , we obtain a set of equations 
which for the case of the central potentials has the form: 

^ + ^{7^(^)-^ A + 1 )} + ^{^ 2 37i^)- / ( / + 1 )}- 
2M k pi I dp k dp k ) 2pjrj 3 I dr j3 J dr j3 ) 

V js )^MM(Pk,r jS ) = I d Pk f dr j3 J2wS M ^ MX/ll ( Pj ,r k3 ), (20) 
where the following notation has been introduced: 

WS M = {Y x (p k ) ® Yfif^Y^ (V j3 + V jk ) {Y x> {p 3 ) ® Mf k ,)} LM ■ (21) 

To progress from (|2"U|) to one-dimensional equations, we apply a modified close coupling 
method, which consists of expanding each component of the wave function ^i(fj 3 ,p k ) over 
the Hamiltonian eigenf unctions of subsystems: 

^3 = ~V^, + V is (53)- ( 22 ) 



Thus, following expansions can be applied: 

^MXi(Pk,r j3 ) = ^J2fn!x M ^) R nhr 3 s), (23) 

P k n 

where functions R l n i(rj 3 ) are defined by the following equation: 

( E » + 2^{4 ( ''^ ) -' (l + 1} } " ^)*"' fe) = °- P4) 

Substituting Eq. ( 1231) into ( 120|) . multiplying by the corresponding functions R l n i(rj 3 ) and in- 
tegrating over r"j 3 drj 3 yields a set of integral-differential equations for the unknown functions 
fnix(Pk)'- 

2M k (E - Ei)r a ( Pk ) + (-fj- ^i^)r a (Pk) = 

v °pk Pk ' 

2M k J2 r dr 3 ,r% f df j3 I dp k -Q^£(pj), (25) 

a JO J J Pj 

where 

Qt' = RUr j3 )wS M Ri ni (r k3 ). (26) 

For brevity one can denote a = nlX (a' = n'l'X), and omit LM because all functions have 
to be the same. The functions f l a (pk) depend on the scalar argument, but this set is still 
not one-dimensional, as formulas in different frames of the Jacobi coordinates: 

Pi = 5'3 ~ PA, *jf3 = -{fikPk + Pj), r jk = -(a-jPj - cr k p k ), (27) 

with the following mass coefficients: 

& = m fc /(l + m fc ), <r k = l-P k , ~? = 1 - Mj (j^k = l,2), (28) 

clearly demonstrate that the modulus of pj depends on two vectors, over which integration on 
the right-hand sides is accomplished: pj = jr^—^kPk- Therefore, to obtain one-dimensional 
integral-differential equations, corresponding to equations (125|) . we will proceed with the 
integration over variables {pj, p k }, rather than {7^3, p k }. The Jacobian of this transformation 
is 7 -3 . Thus, we arrive at a set of one-dimensional integral-differential equations: 

2M k (E - K)P a (p k ) +(-fj~ ^P±)fM = —,Y, (" *Pi&(Pi,Pk)&{Pi)> (29) 

X(J Pk Pk ' 1 a , Jo 

where functions S™ a ,(pj,p k ) are defined as follows: 

Saa'(Pj,Pk) = ZpjPk J dfa I dp k B^(r j3 ) {Y x ( Pk ) ® Yi{f j3 )}* LM (v j3 + V jk ^ 

{Y y (p j )®Y l ,(r k 3)} LM Ri'Ark3) ■ 



One can show (see Appendix, Sect. IV) that fourfold multiple integration in equations 
(130]) leads to a one-dimensional integral and the expression (1301) could be determined for any 
orbital momentum value L: 

S*,(pj,Pk) = 2^ I [(2A + 1)(2A' + l)]hjPk J dwrinwi^fas^Visfas) + 

fy(r iJk ))i&,(7*0 ^^ m ,(0,a; ) 0)(7 A i 7 m ^: m ,y to (^ ) 7r)^ m ,( i / fe ,7r) , (31) 



where D^ m ,(0, u, 0) are Wigner functions, u is the angle between pj and p k , Vj is the angle 
between f k3 and pj, v k is the angle between fj 3 and p k (see the Fig. 1). Finally, we obtain an 
infinite set of coupled integral-differential equations for the unknown functions f^ipi) and 



pap2) m 



/,, n2 d 2 A(A + 1)\ if . ^ (2A + 1)(2A' + 1) f°° , .. . 

/ dusinw-R^^) ( V^r^) + ^'(r ii /)J J R^ / ,(rj 3 )p i /pi 

E £>w(0, w, 0)C A L T m C A ^; m ,y, m (^, vr)^ m ,(^, tt). (32) 

mm' 

For the sake of simplicity a = (n/A) are quantum numbers of a three-body state and L is 



the total angular momentum of the three-body system, g = AnMi/'y 3 , k % n = \j2Mi(E — E%), 
where E l n is the binding energy of the subsystem (i'3), Mi = m 1 (m 2 + 7Ti 3 )/(m 1 + m 2 + m 3 ) 
and M 2 = m2{mi + 777.3)/ (mi + m 2 + 7713) are the reduced masses, 7 = 1 — mj?7v/((mj + 



l)(mj/ + 1)), D„ rm /(0,cj, 0) the Wigner functions, Cj^ m the Clebsh-Gordon coefficients, Yj 



/m 



are the spherical functions, 00 is the angle between the Jacobi coordinates pi and /v, z^ 
is the angle between ry 3 and pi, u^ is the angle between r i3 and /fy. One can show that: 
sinz/j = (p k r kj ) /j sin u, and cos^ = ((3pi + p k cosuj)/{^r kj ). 

C. Boundary conditions, cross sections and numerical implementation 

To find a unique solution to Eqs. (|32|) appropriate boundary conditions depending on 
the specific physical situation need to be considered. First we impose: 

/i°(0)~0. (33) 

Next, for the three-body charge-transfer problems we apply the well known K— matrix for- 
malism. This method has already been applied for solution of three-body problems in the 

10 



framework of the coordinate space Faddeev equations 4l| . For the present scattering prob- 
lem with i + (j3) as the initial state, in the asymptotic region, it takes two solutions to 
Eq.( 152"|) to satisfy the following boundary conditions: 

p^+oo (34) 

fn(Pj) ~ ^/v~/v~ j K ij cos{k < ? ) p j ) , 

where #y are the appropriate coefficients, and V{ [i = 1, 2) is a velocity in channel 2. With 
the following change of variables in Eq. fl32|) : 

f2(P,) = /i ( ?(p,)-sin(A ; f ) p,), (35) 

(i=l, 2) we get two sets of inhomogeneous equations which are solved numerically. The 
coefficients K^ can be obtained from a numerical solution of the FH-type equations. The 
cross sections are given by the following expression: 



47T 
<7, 



V - ,(i)2 



A-r 



K 



zK 



4. M>' + *3 



^ 2 ) 2 ' 



kf 2 (D - l) 2 + (#n + # 2: 

where (z,j = 1,2) refer to the two channels and D = K\\Ki2 — KiiKii- Also, from the 
quantum-mechanical unitarity principle one can derive that the scattering matrix K = 

has the following important feature: 

#21 #22/ 

#12 = #21- (37) 

In this work the relationship ( 1371) is checked for all considered collision energies in both 
antihydrogen cases, i.e. in p + (e + e~) and in p + (fi + fi~), and in the case of the muon 
transfer reactions. 

As stated in Sec lII Al the solution of the Eqs. (J7])-(JE1) involving both components ^1(2) 
required that we apply the expansions (TT9l) and (|23|) over the angle and the distance variables 
respectively. However, to obtain a numerical solution for the set of coupled Eqs. fl32l) we 
only include the -s and -p waves in the expansion (fT0|) and limit n up to 2 in the Eq. 
( )23|) . As a result we arrive at a truncated set of six coupled integral-differential equations, 
since in \&i(2) only Is, 2s and 2p target two-body atomic wave-functions are included. This 
method represents a modified version of the close coupling approximation with six expansion 
functions. The set of truncated integral-differential Eqs. fl32|) is solved by a discretization 

11 



49| and t 



procedure, i.e. on the right side of 

by sums using the trapezoidal rule 

side are discretized using a three-point rule 

equations for the unknown coefficients fa (k) (k — l,N p 



;he equations the integrals over p\ and p 2 are replaced 



re second order partial derivatives on the left 



49] . By this means we obtain a set of linear 



k^ 2 + Dl 



A(A 



>j 



Pu 



./; 



VHi\ - 



M,^^ 



:(i2) 



(2), 



'V-. 



E Z^ w jS aal ) (pu,P2j)f a ,{j) = o 



a'=l j=l 



(38) 



M2 



N s N p 



EE^ 5 ^^'^)^ 1 ^') + 



a=l j=l 



^ + 4 



A'(A' 



p\i 



fS'ii) 



BlXi) 



(39) 



1,AU, N„ is 



<pji 



Here, coefficients Wj are weights of the integration points pu and p 2 i 

the number of quantum states which are taken into account in the expansion (1231) . Next, 

Dfj is the three-point numerical approximation for the second order differential operator: 

Dim = 

A = p i+ i - pi. The vector B^(i) is: B [ J 



(f a (i - l)^_i,j - 2f a (i)5 id + f a (i + i)S i+1>j )/A, where A is a step of the grid 



(21), 



^^ r D 



?(21) 



= M a rf^ Wj SS^(i,j)Bm(k lPj ), and 
in symbolic-operator notations the set of linear Eqs. fl38l) - fl39l) has the following form: 



^ 



^2c/=i Ei=i A. aa i(i, j)f a '(j) = bJ[i). The discretized equations are subsequently solved by 



the Gauss elimination method 



50j | . As can be seen from Eqs. (|38|)-(l39l) the matrix A should 



have a so-called block-structure: there are four main blocks in the matrix: two of them 
related to the differential operators and other two to the integral operators. Each of these 
blocks should have sub-blocks depending on the quantum numbers a = nlX and a' = n'l'X . 
The second order differential operators produce three-diagonal sub-matrixes [47(. However, 
there is no need to keep the whole matrix A in computer's operating (fast) memory. The 
following optimization procedure shows that it would be possible to reduce the memory 
usage by at least four times. Indeed, the numerical equations (138|) - (!39|) can be written in the 
following way: DJ 1 - M^S 12 ] 2 = 0, and - M 2l ~ 3 S 21 f 1 + D 2 f 2 = b. Here, D u D 2 , S 12 
and S 21 are sub-matrixes of A. Now one can determine that: f 1 = (.D 1 )~ 1 M 1 /7 3 S' 12 /' 2 , where 
(Di) _1 is reverse matrix of D\. Thereby one can obtain a reduced set of linear equations 
which are used to perform the calculations: [D 2 - M1M27 -6 S 21 (-D1) -1 S 12 ] f 2 = b [47]. 

To solve the coupled integral-differential equations ( 132]) one needs to first compute the 
angular integrals Eqs. ( 13TT) . They are independent of energy E. Therefore, one needs to 
compute them only once and then store them on a computer's hard drive (or solid state drive) 
to support future computation of other observables, i.e. the charge-transfer cross-sections 
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at different collision energies. The sub-integral expressions in (l31~j) have a very strong and 
complicated dependence on the Jacobi coordinates pi and pi>. To calculate S^ a , .(/?$, pi>) 
at different values of pi and /fy an adaptable algorithm has been applied together with 
the following mathematical substitution: cos a; = (x 2 — f3fpf — p 2 ) / '(20^^) . The angle 
dependent part of the equation can be written as the following one- dimensional integral: 



4tt[(2A+1)(2A' + 1)]3 



ew)= ^ A ^ i y iJJ **2<* 



A 2L + 1 



Pipi+Pi' 



\PiPi-Pi>\ 



-i + 



ru'(x) 



R%l>(r a {x)) 



J2 DL*(0, u(z), 0)C^ m C^;, m ,Y lm (Mx), ir)Y? m ,Mx), n). (40) 



We used a special adaptive FORTRAN subroutine from the work [51| in order to carry out 
the angle integration in ( I40p . This recursive computer pro gram, QUADREC, is a better, 
modified version of the well known program QUANC8 [50[. QUADRE C p rovides a much 



higher quality, stable and more precise integration than does QUANC8 51] . Therefore, our 



results for the three-particle muon transfer reactions presented in Table I are slightly different 



from those of our older work |46|] where we used the less effective adaptive quadrature code 
QUANC8 for numerical computation of the angle integrals ( l4"0l) . The difference between 
these two results ranges from ~9% to ~15%. The expression (1401) differs from zero only in 
a narrow strip, i.e. when p t « pi>. This is because in the considered three-body systems 
the coefficient /3j is approximately equal to one. Figures [3] and H] show the angle integral 
2-dimensional functions (surfaces) (140]) for the (p p~ p + ) system considered herein. For 
example, this might involve a few selected atomic/muonic transitions such as: S[ s .{ s ,(pi, p 2 ) 
and Sl s: 2 S ,(pi,p2). Only the input channel p + (p~p + ) of the reaction (J2J) potential surfaces 
are included. It is seen, that these surfaces have significantly different geometrical shapes 
and numerical values. Therefore, in order to obtain numerically reliable converged results 
it is necessary to adequately distribute a very large number of discretization points (up 
to 2200) between and ~90 atomic/muonic units. More points are taken near the origin 
where the interaction potentials are large and a smaller number of points are needed at 
larger distances. 

III. RESULTS AND CONCLUSIONS 

In this section we report our computational results. Five different three-body Coulomb 
systems have been computed in the framework of a unique quantum-mechanical method, 
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i.e. the FH-type equation formalism ©-((8]). The few-body approach has been presented in 
previous sections. In order to solve the coupled equations ©-flU) we use two different and 
independent sets of target expansion functions (1231) . This is shown in Fig. 2 for the case of 
the (/x + /i~) and (p/i + ) targets. Together with the specific structure of the two coupled FH- 
type equations in the operator form this method allows us to avoid the over-completeness 
problem and the two targets are treated equivalently. The main goal of this work is to 
carry out a reliable quantum-mechanical computation of the formation cross sections and 
corresponding rates of the H and H^ atoms at very low collision energies, i.e. reactions 
([I]) and (J2J). However, as a test of the method and our FORTRAN code we carried out 
calculations of the three-body cross sections and rates of the \i~ transfer reactions from 
d to t: t + (d/i~) — )■ (t/i^) + d, from p to t: t + (p/i~) — > (t/i~) + p, and from p to 
d: d + (p/i~) — > (d/T~) + p. Here, p, d and t are the hydrogen isotopes: proton p= 1 H + , 
deuterium d= 2 H + , and tritium t= 3 H + . The coupled integral-differential Eqs. (132]) have been 
solved numerically for the case of the total angular momentum L = within the two-level 
2x(ls), four-level 2x(ls+2s), and six-level 2x(ls+2s+2p) close coupling approximations in 
Eq. ( 1231) . The sign "2x" indicates that two different sets of expansion functions are applied. 
Next, the following boundary conditions (133|) . (1341) . and (1351) have been used. To compute 
the charge transfer cross sections the expression (1361) has been applied. 

It would be useful to make a comment about the behaviour of a tr (e co u) at very low 
collision energies: e co u ~ 0. From our calculation we found that the muon transfer cross 
sections a tr — > oo as e co u — > 0. However, the muon transfer rates, \ tr , are proportional to the 
product a tr x v c , m , and this trends to a finite value as v c . m . — >■ 0. Here v c . m . = \/2e co u/Mk is a 
relative center-of-mass velocity between the particles in the input channel of the three-body 
reactions, and Mk is the reduced mass. To compute the muon transfer rate the following 
formula is used: 

\ tT = a tr (e coU -> 0)t> c . m .iV , (41) 

where Nq = 4.25 • 10 22 c.m. 3 is the liquid hydrogen density. 

The Coulomb few-body systems mentioned above are of a significant importance in the 
/iCF cycle [2| . In the literature one can find the results of a variety of different calculations of 
these reactions. We compare our results with some of this data. Table I shows cross sections, 
a^, and corresponding thermal rates, Xt r , for all three muon transfer reactions at low collision 



energies together with some theoretical calculations from older papers 361438 
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431 ] and with 



some experimental data from works [52H58| . Our FH-type equation results shown in Table 
I have been computed within the 2x(ls+2s+2p) approximation in the expansion (123]) for 
both Faddeev-type components. Therefore, it was actually used with up to six expansion 
functions. One can see that our a tr and X tr are in fairly good agreement with previous 
calculations and experimental data for all three muonic systems presented in Table I. The 
largest number of results can be found for the first listed reaction, i.e. d+(t/i _ ), which is 
one of the most important three-particle reactions in the /xCF cycle in cold liquid hydrogen. 
We obtained very good agreement with the experimental data and with some theoretical 
results, except in the case of work J37J. A very good agreement is also obtained for the 
other two reactions: for t+(p/i _ ) and d+(p/i _ ). These results show that the few-body 
method of the two coupled FH-type equation ([7|)-(JH]) works extremely well together with 
the close coupling expansion method (123]) . Only three therms in the expansion (123]) . i.e. 
2x(ls+2s+2p) approximation, can provide such good agreement with the experiments for 
all three muonic transfer reactions. As we already mentioned the sign "2x" means that 
the three term expansion is used within two different expansion functions sets. Figs. 5 
and 6 show our cross sections for t+(p/i~) and d+(p/i~) collisions. Again, these results are 
obtained within the different close coupling expansion approximations: 2xls, 2x(ls+ 2s), 
and 2x(ls+2s+2p). In the first case only two expansion functions are used, in the second 
case only four, and in the third case six expansion functions are used. One can see how 
significantly the 2p-state contributes to the total muon transfer cross section when decreasing 
the collision energy, especially this is seen in the case of d+(p/i~). A comparison with the 
experimental data and other theoretical results in Table I for all three muonic transfer 
reactions demonstrates that the FH-type equations and 2x(ls+2s+2p) approximation are 
able to provide reliable results for three-body charge transfer reactions at low energies. 
It is also important to mention here, that the pure quantum-mechanical behaviour of the 
transfer cross section at low energies, specifically a tr {ScoU ~ 0) ~ l/v c , m ,, has been obtained 
in this calculation. This allowed us to compute the muon transfer rates ( fld]) . i.e. Xt r {T — > 
0) ^sconst, and compare these results with the experiments. 

Next, the three-body reaction of the atomic H formation, i.e. reaction (CQ) is considered. 
We are primarily interested in low energy collisions. Because this is not a muonic system 
one needs to switch from muonic to atomic units. Also, in the computer program one needs 
to change the masses of the particles. Fig. 7 shows our results for the reaction ([1]) total 
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cross section in the framework of different close-coupling approximations in the Eq. (j23|) . 
One can see, that the contribution of the 2p-states in each target become larger while the 
collision energy becomes smaller. As in the muonic transfer reactions we found that the cross 
section of the H formation ctjj — >■ oo as e co u — > 0, i.e. a^v c , m . ~const. This fact allows us to 
compute the low energy rate of the H production. For example, one can follow the logic of 
work |3] and estimate the H production rate by using the following formula: Rg = a^N-plI. 
According to [3] Np is the density of slow antiprotons, / is the number of p traversing the 
interaction region each second, and I is the linear dimension of the interaction region. In 3[ 
the last parameter was taken as I — 1 cm. The product of // has the unit of velocity, thus 
it should be possible to represent the rate R^ as well as the expression ( HIT) , i.e.: 

i?H = ^H Vc.rn.Np, (42) 

where v c , m , is the cm. velocity between p and the positronium atom Ps. Our results for 
H are shown in Table II together with the results for the H^ formation reaction fl2]). At 
low energies R^ starts taking a constant value as was the case in our previous calculation 
of the /i~ transfer reactions. It would be interesting to estimate the Np parameter. For 
example, if we accept the recent data from the Evaporative Cooling (EC) experiment at 



20J: n EC -4000 very cold p at T EC » 9 K, 



the Antiproton Decelerator (AD) at CERN 
we would need to place this quantity of antiprotons in a limited space with a volume V^ c ■ 
The temperature Tec corresponds to the low energy collisions considered in this work: 
~ 10~ 4 eV. In order to obtain the rate R^ > 1, one would need the following volume: 
V^ c » 4000cJh v c . m . = 4000 x 6.68 x 10" 9 « 27 x 10" 6 cm 3 , where the value for a^ v c . m . is 
taken from Table II. If we suppose that the interaction region between p and the Ps atoms has 
a cylindrical shape with the length Iq = 1 cm [3[, its radius r$ should be: tq ~ 8.6 x 10~ 3 cm. 
As one can see r has a very small value, although it seems to us that it still would be possible 
to adopt this value in some experiments. However, a different situation arises if we adopt the 



results of a newer experiment on Adiabatic Cooling (AC) of antiprotons 2l|. The authors 
of this work obtained 3 x 10 6 cold antiprotons at temperature 3.5 K! In this case one would 
need the following volume: V^ c »3x lOVjj v c . m . = 3 x 10 6 x 6.68 x 10~ 9 ps 20 x 10~ 3 cm 3 . 
It is easy to compute that in this case the radius is tq ~ 0.1 cm. Finally, Fig. 8 shows our 
results for the H M formation cross section. It is clear that in the process ([2]) the contribution 
of the 2p-states from each target is becoming even more significant while the collision energy 
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becomes smaller. Additionally, for the process ([2]) we also compute the numerical value of 
the quantity: oyf (e co u — > 0)t> c . m . ~ const. Table II includes our data for this important 
parameter together with the H^ formation total cross section. All these results are obtained 
in the framework of the 2x(ls+2s+2p) close coupling approximation. These data can be 
useful in future developments of low energy collision experiments with participation of cold 
antiprotons and true muonium atoms. Next, because of the complexity of the few-body 
method, in this work only the total orbital momentum L = has been taken into account. 
It was adequate in the case of slow and ultraslow collisions discussed above. However, to take 
into account the important contribution of higher L's at higher collision energies it would be 



possible to use Takayanig's Modified Wave Number Approximation (MWNA) method 59 ]. 



In the recent work [60[ the MWNA method has been successfully applied to a few-body 
charge transfer reaction. In conclusion, it is feasible to expect that the FH-type equation 
formalism ©-flU) could also be an effective tool for computation of the quite intriguing three- 
body reaction ([3]). This is another process of the H M atom production with a heavy charge 



transfer from one center to another 
rates of both processes (J2j) and (J3J 
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3l| . It would be interesting to compare the reaction 



hat 



13| 



221 ] . An additional point to emphasize would be t 
in some sense the reaction ([3]) and the few-body protonium (Pn) formation reaction 9- 
show close similarities. Therefore, it would be good to try to apply the FH-type equation 
method to the Pn formation problem too. Also, it seems quite possible to expand in some 
way or another the FH-type few-body equation approach and the modified close-coupling 
expansion method, Eqs. ( 1191) and ( |23l) . to very important but challenging low-energy four- 
particle rearrangement scattering collisions with the pure Coulomb interaction between the 
particles, such as H+H— > (pp + ) n i + (e + e~) n 'i/ or, for example, H^+H^ — > (pp + ) n i + (/U + / i ~)n'i'- 

IV. APPENDIX 

The details of the derivation of the angular integrals S™ a ,(pj, pk) (T3TT) are explained below 
in this section. The configuration triangle, A(123), is determined by the Jacobi vectors 
(r*j3,pfc) and should be considered in an arbitrary coordinate system OXYZ. In this initial 
system the angle variables of the three-body Jacobi vectors {rj3,Pk} have the following 
values: fj 3 = (9j,<f)j), pk = (Qk,$k), j 7^ k = 1,2. Let us adopt a new coordinate system 
O'X'Y'Z' in which the axis O'Z' is directed over the vector pk, A(123) belongs to the plain 
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O'X'Z' and the vertex k = 1 of A (123) coincides with the origin O' of the new OX'Y'Z' . 
Fig. 9 shows the specific configuration of A (123) and the new adopted O'X'Y'Z' system. 
One can see, that the new angle variables of the Jacobi vectors in the O'X'Y'Z' system 
have now the following values: f'j 3 = (ffc,7r), p' k = (0,0), r' k3 = (r] k ,ii), p'j = (oj,it), here 
k — 1 and j = 2. The spatial rotational transformation from OXYZ to O'X'Y'Z' has been 



done with the use of the following Euler angles (Q k ,Q k ,e) 48j. Taking into account the 



transformation rule for the bipolar harmonics between new and old coordinate systems, one 



can write down the following relationships 48j: 



{yx{pk)®Y l {? J s)} m LU = J2( D M m (®k,e k ,e)y {Y x (p' k ) ® i^K™ (43) 

in 

{Y xl {p 3 )®Mf k ,)} LM = J2 D M m >(®k,e k ,e) {Y A ,(#) ®^(^ 3 )} w , (44) 



where D\ lm {Q k) 0fc, e) are the Wigner functions [48]. The fourfold multiple angular inte- 
gration j dpj j df>k in Eq. (I3"0"|) can be written in the new variables and be symbolically 
represented as j dpj j dp k = f dui sin uj J Q n de f Q d<& k Jq sin Q k dQ k . Next, taking into ac- 



48 



count the normalizing condition for the Wigner functions 

r'I'K r'I'K r>n 8*77-2 

/ de d$ k sme k de k (D^ m (^ k ,e k ,e)YD^ m ,(^ k ,e k ,e)) = ^-^5 mm i (45) 
one can obtain the following intermediate expression: 

(V ja + V jk ) {Yxift) ® Y v (f' k3 )} Lml RiAns) ■ (46) 

Now, let us make the next transformation of A(123) in which the vertex j = 2 of A(123) 
coincides with the centre O' of the O'X'Y'Z' and O'XYZ, however the axis O'Z" is directed 
along pj and A (123) belongs to the plain O'X" Z" . This transformation, which converts the 
coordinate frame O'X'Y'Z' into 0'X"Y"Z" is characterized by the following Euler angles 
(0, u, 0). Therefore the vectors (r fe3 , pj) have the following new variables: f k3 = (fj, n), p" = 
(0, 0). As a result of this rotation one can write down the following relationship: 

{YAp'j)®Yi'(r' k3 )} Lm = J2 D Mm^,^0) {YAp")®Yv{r'U)} Lm , (47) 

in 1 

and obtain the following result: 

S % L{P 3 ,Pk) = 2p jPk J2 ^7 J duJ sin <"J&(r is ) i Y ^ °) ® Y ^)Y Lm Ok + V jk ) 

^w(0, w, 0) {F v (0, 0) ® Mf£ 3 )} w i&,(r M ).(48) 
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Now by taking into account that YJ m (0,0) = 8 m $y/(2l + l)/47r 48], the bipolar harmonics 
in (SHI) are: 



{F A (0, 0) ® Y t (u k , n)Y Lm = yj^lcfZXim^, *), (49) 

/oTi _i_ i 
{Fv(0,0)®^(^^)}w = y— £-(%$,„,*,„,&,*), (50) 

with the use of these relationships we finally get the convenient for numerical computations 
Eq. (|3H). 
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TABLE I. Cross sections at T and rates \t r , Eq. (|4ip . for fj,~ transfer reactions from a light hydrogen isotope to a heavier hydrogen isotope 
at low collision energies together with other theoretical results and experimental data. The result for unitarity ratio K21/K12, Eq. (j37j) . are 
also presented for t+(d^)i s . 



Energy, eV Method t+(d/i)i s -> (t^)i 3 +d 

a tr /W- 20 , cm 2 A tr /10 8 , s" 1 



t+(PM)ls -» (t/i)ls+P 



d+(PA0ls -> (d/i)i s +p 



K 2 i/K 12 



vtr/W- 20 , cm 2 



A tr ./10 8 , 



<Ttr/10- 20 , 



A tr /10 8 , 



to 

10 



0.001 



FH-type Eqs.: 
[36] 

[37] 

[38] 

[39] 

Experiments: 



15.4 
15.8 
21.5 
18.0 
14.2 



2.6 
2.7 
3.5 

2.8 

2.8±0.5[58J 
2.8±0.3[57J 
2.9±0.4[56J 



0.99 



315.0 


65.1 


384.4 


80.0 


265.0 


55.0 



663.4 


146.0 


828.7 


170.0 


650.0 


140.0 





58.6±10[55] 




84±13[54J 

143±13[53] 

95±34[5J] 


99.1 


64.7 


208.1 


144.8 


128.0 




283.7 




60.0 




140.0 





0.01 



FH-type Eqs. 
[36] 

[37] 
[38] 
[39] 



4.84 
5.64 

1.8 
5.0 
4.44 



2.6 



0.99 



0.04 



FH-type Eqs. 
[36] 

[37] 
[38] 



2.37 
2.94 

3.1 

2.5 



2.5 



0.99 



49.1 
63.6 
40.0 



64.2 



103.1 
140.7 
91.0 



143.4 



0.1 



FH-type Eqs. 
[36] 

[39] 



1.42 

2.0 

1.35 



2.1 



0.99 



30.7 
39.9 



63.5 



64.6 
87.4 



142.1 



TABLE II. The total cross sections crjj and <tjj for the reactions ([I]) and ([2]) respectively. The 
product of these cross sections and the corresponding center-of-mass velocities u c .m. between p and 
Ps=(e + e~), i.e. cr^v c , m , and between p and the true muonium atom Ps M = (fj, + fi~), i.e. ctjj v c . m , 
are presented. 



p + (e + e ) ls ->H + e~ 



p + <yv )i s -> H M + ^' 



£, eV 



ojj-, cm 



ffffcm., cm 3 /s 



(7 H cm 



cgWcm.! cm 3 /s 



1.0e-06 


0.16e-12 


0.67e-08 


1.0e-05 


0.50e-13 


0.67e-08 


1.0e-04 


0.16e-13 


0.67e-08 


1.0e-03 


0.50e-14 


0.66e-08 


1.0e-02 


0.15e-14 


0.63e-08 


5.0e-02 


0.60e-15 


0.56e-08 


1.0e-01 


0.42e-15 


0.55e-08 


5.0e-01 






1.0e-00 







0.18e-16 
0.58e-17 
0.18e-17 
0.82e-18 
0.58e-18 
0.27e-18 
0.23e-18 



0.60e-12 
0.60e-12 
0.59e-12 
0.59e-12 
0.59e-12 
0.62e-12 
0.73e-12 
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FIG. 1. The configurational triangle of a three-charged-particle system (123). (p/u~/x + ) is presented 
together with the Jacobi coordinates, the inter-particle vectors, and the angles between the vectors 
and coordinates. O is the center of mass of the few-body system system, 0\ and O2 are the center 
of masses of the targets \i~ \x + and p// + respectively. 
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FIG. 2. Two asymptotic spacial configurations of the 3-body system (123), or more specifically 
(p, fjT , fi + ) which is considered in this work. The few-body Jacobi coordinates (pi,fjk), where 
i/j/l:= 1, 2, 3 are also shown together with the 3-body wave function components ^i and ^2' 
^ = \&i + ^2 is the total wave function of the 3-body system. 
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FIG. 3. The angular integral S™ a >(pi,Pi'), Eq. PU|) . in the input channel p + (fi + fj, ) when a =ls 
and a' =ls. 



26 










FIG. 4. The angular integral S™ a >(pi,Pi'), Eq. PU|) . in the input channel p + (fi + fj, ) when a =ls 
and a! =2s. 
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FIG. 5. Low energy cross sections for a muon transfer reaction in the t+(p// _ )i s collision. Results 
are shown for the two-level 2xls, four-level 2x(ls+2s), and six-level 2x(ls+2s-|-2p) close-coupling 
approximations. 
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FIG. 6. Low energy cross sections for a muon transfer reaction in the d+(p/j,~)i s collision. Results 
are shown for the two-level 2xls, four-level 2x(ls+2s), and six- level 2x(ls+2s+2p) close-coupling 
approximations. 
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FIG. 7. Low energy cross sections for an antihydrogen atom formation reaction ([I]). Results are 
shown for the two-level 2xls, four- level 2x(ls+2s), and six-level 2x(ls+2s-|-2p) close-coupling 
approximations. 
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FIG. 8. Low energy cross sections for a muonic antihydrogen formation reaction ([2]). Results 
are shown for the two-level 2xls, four-level 2x(ls+2s), and six- level 2x(ls+2s-|-2p) close-coupling 
approximations . 
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FIG. 9. The configurational A(123) in the case of the p + (fj,~ fi + ) collision is shown together with 
the new O'X'Y'Z' cartesian coordinate system after the rotational-translational transformation 
from the initial OXYZ system (see Appendix, Sect. IV). OXYZ is not shown here. 
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